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A gauge invariant electromagnetic two-point function, crucial to the inves- 
i-£h ' tigations of the violation of isospin symmetry, is derived for heavy-light quark 

systems. Thus QED can be consistently introduced into the QCD Sum Rule 



Oh 
D 



Q-i. method, which was not previously possible. 
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^ . The origin of mass differences in isospin multiplets has long been of great interest 

in nuclear and particle physics as a source of information about symmetry violations. 
Hadronic isospin violations are particularly important in that they arise from non- 
perturbative Quantum Chromodynamics (QCD) as well as quark mass differences 
(see Ref.[0] for a review of the early work in this area), and of course electromag- 
netic effects. Among the first applications of the method of QCD sum rules was 
the study of isospin violations in the p — oj system 0, in which it was recognized 
that the isospin splitting of the light-quark condensates can produce effects as large 
as the current-quark mass splittings and electromagnetic effects. Recently the QCD 
sum rule method has been used to study the neutron-proton mass difference [§, |J, 
the octet baryon mass splittings and the mass differences in the charmed meson 
systems (the D and D* scalar and vector mesons) @. 

In comparison with hadronic quark models, the QCD sum rule method for calcu- 
lating isospin splittings has the advantage that one can directly use QED field theory 
rather than rely on models to estimate Coulomb corrections. This, however, has not 
been done. In our earlier attempt to calculate the standard two-loop QED contribu- 
tions to the isospin mass splittings in heavy-light quark systems [Fig. 1] using the 
sum rule method we found |7| that for the charged mesons, involving charged currents, 
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the calculation is not gauge invariant. The objective of the present letter is to devel- 
ope a gauge invariant theory for QED phenomena within the QCD sum rule method. 
Although our derivation is for heavy-light quark mesons, the method is quite general 
and can be used for baryons as well as mesons. 

The basic approach of the QCD sum rule in heavy light quark systems is to study 
the two-point function in the Wilson operator product expansion (OPE), defined by 



KM 2 



d 4 xe iqx < T(J„(x)J v (0)) > 



(1) 



for the heavy-light quark current 



(2) 



The local operators, {O n }, consisting of quark and gluon fields and the Wilson coef- 
ficients functions, |/^(q 2 )}, have been extensively discussed in the literature in the 
studies of the masses and their decay constants for heavy-light quark systems. 

In order to study the violations of the isospin symmetry, the electromagnetic effects 
should also be written in the framework of the operator product expansion. The 
leading electromagnetic effects in the this approach are the two-point functions from 
a two- loop perturbative contribution, whose Feynman diagrams are shown in Fig. 1. 
For the charge neutral current, one could simply obtain the two-point functions by 
changing the gluons in QCD to the photons in QED in Fig. 1, since the two-point 
functions has been calculated in QCD|§. Their imaginary parts are 
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for pseudoscalar and vector currents, where x = 



9 /5 
f(x) = - + 2l(x) + ln(x) ln(l - x) + I - - x - — 
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X 
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— — x ) ln(l — x) (5) 



and l(x) = — J(fln(l — y)-^ is the Spencer function. Obviously, this result is gauge 
invariant. 

The calculation of the two-point functions for a charged current is much more 
complicated. Since the Feynman diagrams in Fig. 1 are weighted by the products of 
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different charges, the corresponding two point functions are not gauge invariant. The 
physical origin of this problem is that under the gauge transformation 



q(x) -> e ie " A{x) q(x) 
Q(x) -> e ie « A(x) Q(x) 
A^x) -»• Ap(x) -d»A(x), (6) 

we have 

J Ax) - = e iA ^^ +e ^g(x) 7M Q(x) = e ^K e Q+%) J^z) (7) 

(notice e g = — e 5 ), therefore, the current becomes gauge dependent with e q + 
e Q — e r- Of course, gauge invariance is a fundamental property of electromagnetic 
interactions, and any observable obtained using the current must be gauge invariant. 

In order to obtain a solution to the problem, let us start with a gauge invariant 
form of the two-point function, which is obtained by inserting the link operator into 
the form given in Eq. [I], with the definition: 

IV(g 2 ) = i J d'xe^-v) < T(J,(x)e lQ °^y Aa(y)dya l(y)) >, 

= iJd 4 xe i ^<T(Jl(x)Ji(y))> (8) 

where Q op is the charge operator and 

Jj(a;) = q(x)e iQ ° p fo A ^ d y a ^Q( x ). (9) 

The gauge transformation of the electromagnetic field A a (x) cancels those of quark 
fields q(x) and Q(x) so that two-point function is gauge invariant. This is QED analog 
to the QCD treatment introduced for the pion wave function. 

Expanding to order a e one finds that there are two currents at the two loop level: 

4(x) = j°(x) + j;(x). (io) 

The current J®(x) is given by insertion of photon vertices in the quark lines of the 
current given by Eq. Q and the evaluation of the resulting two-point functions is 
carried out by the standard two- loop diagrams shown in Fig. 1. The charges involved 
are eQ and e q for these processes. The second current of Eq. [TOj , 

J^x) = ie T q(x)^ J A a (y)dy a Q(x), (11) 

corresponds to an additional vertex function with the charge e^; and in the evaluation 
of the two-point function gives rise to the additional diagrams shown in Fig. 2, which 
we now discuss in detail. The resulting gauge invariant two point function is 

= K (? 2 ) + nU<z 2 ), (12) 



3 



where 

KM 2 ) = * J d 4 xe^ < T (Jj(x)^(0)) > (13) 
at the two-loop level is given by the Feynman diagrams in Fig. 1, and 

= %J d*xe^r {j°(x)r u (o) + j;(x)j° u (o) + r,(x)r u (o)} (u) 

generates additional Feynman diagrams shown in Fig. 2. These additional terms 
introduced by arise from the additional vertex function corresponding to J^. 
The fact that such additional vertex functions must be present for charged currents 



for gauge invariance has been observed previously. See, e.g., Ref.jlO 



Since the masses for the heavy and light quarks are not the same, it is more 
convenient to calculate the two-point function in momentum space. The Fourier 
transformation of the operator Jq A^^dy 11 gives 



F{q 2 ) = J tfxji* £ A^dy** 



where g M is the momentum carried by photon field A^(q 2 ), therefore the current J^{q 2 ) 
in the momentum space is given by 



J e M 2 ) = -e T q(q ~ h) % || _ ^ ^{h - k 2 )Q{k 2 ). (16) 



The sum of the Feynman diagrams in Fig. 2 has a simple form 

e T {e q + e Q ) f d D k 1 d D k 2 Tr[^^^ 2 ly{h 



TP (n 2 ) 6T[6q + 6q) f 

A VW ) 28yr 4 J 



n D {h - k 2 )%k% - M 2 )(h - q) 2 
= e T (e q + e Q )I^(q 2 ). (17) 

To show that the addition of H e (q 2 ) indeed makes the total two-point function, 
fl^g 2 ), gauge invariant, we separate fl^g 2 ) into the gauge dependent and indepen- 
dent parts; 



KM 2 ) = e 2 ^ u (q 2 ) + e 2 Q I^q 2 )-e q e Q I^(q 2 

= -e q e Q (ilM 2 ) + i?M 2 ) + ifAq 2 ) 

+e T e q I%(q 2 ) + e T e Q I^(q 2 ), (18) 

where I^iq 2 ) and I® u (q 2 ) represent the the self-energy diagram for the light quark q(x) 
and heavy quark Q(x) respectively, and I^{q 2 ) corresponds to the photon exchange 
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between the light and heavy quarks in Fig. 1. The term with I^iq 2 ) + lf* u (q 2 ) + 
Ifa!(q 2 ) in Eq. [18] is of course gauge invariant, and the analytical expressions of its 
imaginary part can be obtained from Eqs. |3] and The last two terms in Eq. [IB] are 
gauge dependent; if we substitute the gauge dependent part of the photon propagator 
( D ^) g d = (kl ~tkl-^ k2)v into the two-loop integrals I^(q 2 ) and 1% (q 2 ), the gauge 
dependent part of the integrals is 

fetf)),, = = - Uo 2 ) (W) 

where I^ y {q 2 ) is given in Eq. [T7|. By substituting the gauge dependent photon prop- 
agator into n^,(g 2 ) one can show that the gauge dependent part of the two- loop 
integral for H e (q 2 ) is identical to the expression in Eq. [T7]. Therefore, the gauge 
dependent two-loop integrals in 11° (q 2 ) and Tl e (q 2 ) exactly cancel each other, and 
we have a total gauge invariant results. This provides an important check at the 
two- loop level that the two point function (q 2 ) for the current Jt(x) in Eq. |l(] is 
indeed gauge invariant. 

To evaluate the two point function Il^ u (q 2 ), we rewrite Eq. [1^ as 

iW) = nfAq 2 ) + n^(g 2 ) + ng, (g 2 ), (20) 

where 

(q 2 ) = ~e q e Q (l« u (q 2 ) + I%(q 2 ) + I^(q 2 )) (21) 
whose analytical expressions can be obtained from Eqs. ^| and |], 

K^ 2 ) = ^e q {llM 2 ) + I,M 2 )) (22) 



and 

For Il^(g 2 ), we have 



K^) = ^e Q {l%{q 2 ) + lM))- (23) 



VW ) + V19 ) ~ 2 s % a J (A .2 _ M2)(A . 2 _ q y > ^) 



where 



W - / S Tr^i^ ^ - '^'^H (25) 



7T 



D/2 _ fc 2 )2( fcl _ ? )2 ^ _ fe 



is an one-loop wavefunction renormalization equivalently in Landau gauge. It has 
been shown that vanishes [O] for a zero mass particle in dimensional regular- 

ization. This leads to 

K%{q 2 ) = e T e Q {I%{q 2 ) - I^q 2 )) (26) 
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where (q 2 ) is proportional to the mass of the heavy quarks. This shows that the 



divergence induced by the wavefunction renormalizations does not exist in Eq. |26], 
which implies that the Ward identity is restored in this approach. 

After including the mass renormalization, the two- loop integral U® u (q 2 ) in the 
dimensional regularization is 



47T 



where 



3a P 



-M 



d 



47T6 dM 



W) + 



jO I „2\ 



fiu 



(q 2 



q 2 e^Y f d D k Tr {^{jt + M)j„{jt- 



16vr 2 

is the one-loop integral, D = 4 — 2e 
and the running mass M(/i) is related to the pole mass M by 

,2 

M(/z) =M\l 



(27) 
(28) 



n D/2 (k 2 - M 2 )(k - q) 2 
-> + , is the number of spacetime dimensions, 
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(29) 



The evaluation of Eq. [27] is performed in the modified minimal subtraction scheme 
(MS). Here we only present the imaginary parts of Ti^ v (q 2 ) for the pseudo-scalar and 
the vector states; 



3a e e Q e T M 2 

167T 2 
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(31) 



for pseudo-scalar and vector currents respectively, in which the running mass M(/i) is 
replaced by the pole mass M. In the limit of M — > 0, U®(q 2 ) vanishes for both pseudo- 
scalar and vector currents. Thus, the two point functions are dominated by the term 
proportional to e 9 eg in the small M limit, this is exactly what one would expect 
from the phenomenological model. However, how large effects the contributions from 
n 1 (q |2 ) are in the heavy quark limit remains to be studied. 

In summary, we have derived a gauge-independent form for QED corrections to 
QCD sum rules by identifying processes with additional vertices which must be in- 
cluded for charged currents. The result here presents an important step toward a 
consistent treatment of the violations of isospin symmetry in the QCD sum rule 
framework. The applications of this result to isospin splittings of heavy-light quark 
systems as well as the Kaon systems are in progress, and will be given elsewhere. 

This work is supported by National Science Foundation grant PHY-9023586. 
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Figure Captions 

1. The two loop perturbative corrections for charged neutral currents 

2. Additional diagrams generated by Jt{x) for charged heavy-light quark sys- 
tems. See text. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409377vl 



